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We give the exact solution of orbit dependent nuclear pairing problem between two non-degenerate 
energy levels using the Bethe ansatz technique. Our solution reduces to previously solved cases in the 
appropriate limits including Richardson's treatment of reduced pairing in terms of rational Gaudin 
algebra operators. 
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Pair correlations are manifest in a remarkable range 
of quantum many-body systems. Originally the idea of 
pairing and the methods were developed in the context 
of superconductivity in macroscopic systems by Bardeen, 
Cooper and Schrieffer [l| and by Bogoliubov 0. In nu- 
clear physics, it has been long known that the indepen- 
dent particle picture must be improved with the addition 
of a nondiagonal two-particle force as evidenced by the 
absence of single particle excitations at low energies J3| 
and soon the idea of the pairing was carried over [J, [J. 
But methods borrowed from superconductivity in infi- 
nite systems were inconvenient for the finite nucleus be- 
cause the former is based on wave functions with indefi- 
nite number of particles. Although it is known that the 
pairing effects play a major role in determining nearly 
all nuclear properties including the excitation spectra, 
the transition probabilities and the equilibrium shape, 
an exact number-conserving solution to nuclear pairing 
problem is still lacking except in three limits. The limit 
in which single particle energy levels are degenerate and 
all pairing strengths are the same was solved by Ker- 
man using the quasi-spin algebra @. Later Richardson 
solved the limit in which the single particle energy levels 
are nondegenerate but the pairing strengths are the same 
Q. Finally, the limit in which the single particle levels 
have different pairing strengths but are all degenerate is 
solved by Pan et al [1| and by Balantekin et al [3, ■ On 
the other hand, in many cases it is the interplay between 
the one-body and the two-body effects which determine 
the equilibrium properties of the nucleus and a full solu- 
tion of the nuclear pairing problem is highly desirable as 
a realistic model. The purpose of this Letter is to present 
an exact solution of the pairing problem away from the 
above mentioned limits, i.e., with orbit dependent pair- 
ing strengths and non-degenerate single particle energies, 
in the presence of two nuclear levels. 

The problem is described by the Hamiltonian 



spin operators S^''^ are given by 



\G\Y.c*cyS+S-,. (1) 
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Here c, 's are the occupation amplitudes and the quasi- 
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S^j = ^ - (^a] ^Qj „i + a] -mOj _„i — 1^ . (2) 
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The interaction described by the Hamiltonian ([T]) is 
known as separable pairing, since the interaction strength 
is the multiplication of two numbers (i.e. CjCj'). ^ 

The limit of the pairing problem described by ([T]) 
in which all single particle energy levels are degenerate 
(leading to the first term being a constant for a given 
number of pairs) and all Cj 's are the same can be treated 
using the quasi-spin algebra generated by S^'^ 0]. 

The case when all Cj 's are still the same, but single par- 
ticle energies are non-degenerate was treated by Richard- 
son 01 who obtained energy eigenvalues and eigenstates 
in terms of the rational Gaudin algebra operators 



J+(x) = 



(3) 



dependent on solutions of certain Bethe ansatz equations. 
Here Sj = ej/|G| are the scaled single particle energies. 
More recently, the authors of Ref. ;&] solved the case 
when the single particle energies are degenerate, occupa- 
tion amplitudes of individual orbits (cj's) are different, 
but the shell is at most half full, using the step operators 



-SI 



(4) 



This solution was generalized to the case when the shell 
is more than half full in Refs. @ and (Toj . 



yamac@physics.wisc.edu| 



^ Certain nonseparable pairing models are also known to be inte- 
grable such as the rational, trigonometric and hyperbolic Gaudin 
magnet Hamiltonians [Til . |12|| . These three sets of Hamiltoni- 
ans are mutually commuting and simultaneously diagonalizable. 
They can be combined in various ways to built other integrable 
pairing models. The separable Hamiltonian in ((TJ, however, does 
not commute with any of the Gaudin magnet Hamiltonians and 
hence do not belong to the class of pairing problems which can 
be approached in this manner except when all Cj 's are the same. 
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To give an exact solution for the most general problem 
with two orbits we use the step operators 



J+ix) = 



3 ■' 



(5) 



which approach to those in given ([3]) and ([¥]) in the ap- 
propriate limits with a rescahng of the variable x. 
To present our solution, we write ([1]) as 



^c*c,.^+5-+^e,l]„ (6) 



where we used ([2]) and divided by \G\ to work with scaled 
single particle energies Sj = ej/|G|. In the last term of 



Hamiltonian 



17 o denotes the maximum number of 



pairs that can occupy level j. Although this term is a 
constant, we keep it in order to guarantee that the energy 
of the empty shell is zero. The empty shell is represented 
by the lowest weight state of all the quasispin operators 
defined in ^ and we denote it by |0). 

To find the energy eigenstates with one pair, we form 
the Bethe ansatz state J"'"(a;)|0) where J^{x) is given 
by (O. Using the quasi-spin algebra generated by the 
operators introduced in one can show that the state 
J+(x)|0) is an eigenstate of Hamiltonian ([6]) with the 
energy 



13-x 

if X is chosen so that the Bethe ansatz equation 
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is satisfied. Here a and /3 are given by 



a = 2- 



This result generalizes to the state 



J+{x^)J+{x2)...J+{xNm, 



(7) 



(8) 



(9) 
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with N pairs where all the variables Xk are different from 
one another. Here we assume that the shell is at most 
half full, i.e., N < Nmax/^ where N„iax is the maximum 
number of pairs that can occupy the shell. As will be 
described below, if the shell is more than half full, it is 
easier to work with the hole pairs instead of the particle 
pairs. One can show that the state ifTO]) is an eigenstate 
if the parameters Xk satisfy the Bethe ansatz equations 
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If these equations are satisfied, then the energy of the 
state is given by 



E 



N 



N 
n—1 



(12) 



In general, equations of Bethe ansatz have more than one 
solutions. Each set of numbers {xi,X2, . . . ,xn) satisfying 
the Bethe ansatz equations gives us an eigenstate in the 
form of (jlOp . Solutions may involve complex numbers 
but since they always come in complex conjugate pairs, 
the energy given in ()12|) is always real. 

It is easy to see that the two particular cases, namely 
the orbit independent (i.e., reduced) coupling case solved 
by Richardson and the degenerate single particle energy 
case solved by Pan et al and Balantekin et al can be 
obtained from HIl) and P^ . 

Richardson Limit (/3 —>■ oo): In the limit where the 
occupation amplitudes are all equal to each other (i.e. 
Cj — c), the step operators of (O approach to rational 
Gaudin algebra operators given in ([3]) with a rescaling 
of the variable x a.s x' — |cpx. Consequently, the eigen- 
states (|10|) approach to those proposed by Richardson 
[3]. Since we have /3 ^ oo and a//? — > |cp in this limit, 
Bethe ansatz equations pT|) and the energy ^2]) become 



N 

= l(#fc) " 
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(13) 



(14) 



These equations are those obtained by Richardson in 0] . 

Degenerate Limit{(3 — > 0): In the limit where single 
particle energies are degenerate (i.e. Ej = e), the step 
operators given in ([5]) approach to those in ^ with a 
rescaling of the variable x as x" = x/2e. Therefore the 
eigenstates (ITU)) approach to those found in Refs. [1, H, 
|lO[. Since /3 ^ in this limit, one should distinguish 
between the case in which all variables Xk are different 
from zero and the case in which one of the variables is 
zero. Consequently, we obtain two types of eigenstates 
in this limit. 

From pT|) . it can be easily seen that if all Xk are dif- 
ferent from zero, then the Bethe ansatz of pT|) becomes 



N 



5^ „// 



(15) 
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in the degenerate limit. Corresponding energy given in 
approaches to a constant value given by^ 



En = 2A^e. 



(16) 



Some solutions of the equations of Bethe ansatz pT|) 
involve one variable, say xi, which approaches to zero as 



^ These states are called "zero energy states" in Refs. |^ [lo| 
since the degenerate single particle energy £ is taken to be zero 
in which case the energy JTSJ also vanishes. 
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FIG. 1: Energies for two particle pairs with ji — 3/2 and 
j2 = 5/2. Here i9 describes the occupation amplitudes as 



described in the text and A = 
between single particle levels. 



is the separation 



Equations (fT5llT8| are those found in Refs. 

If the shell is more than half full, it is easier to work 
with hole pairs instead of particle pairs. Fully occupied 
shell is described by the highest weight state of all the 
quasispin operators given in ([2]). We denote this state by 
0). The energy of the fully occupied shell is 



E 



(19) 



Here N^ax denotes the maximum occupancy number of 
the shell. We can create hole pairs (or, equivalently, an- 
nihilate particle pairs) by acting on the fully occupied 
shell with the operators 



J-(y) = E 



2ew 



2 5- 



(20) 



One can show that the state J^{y)\Q) which has one hole 
pairs {Nmax — 1 particle pairs), is an eigenstate of the 
Hamiltonian with the energy 



E 



-i=E_ 



y 



if y obeys the Bethe ansatz equation 
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Similarly, for N hole pairs {N, 
write the Bethe ansatz state 



2£,-|c,fy P-y- 

— N particle pairs) , we 



J-(yi)J-(2/2)...J-(y^)|0). 



(23) 



Here N < Nmax/ "2 and all the variables yk are differ- 
ent from each other. This state is an eigenstate of the 
Hamiltonian ^ with the energy 



N 



E 



v = E, 



_^yn 



(24) 



if the variables yk obey the equations of Bethe ansatz 
given by 



the single particle levels become degenerate. Note that 
there can be at most one such variable since all Xk must 
be different from one another. In this case, the limit 
of the Bethe ansatz equations pT|) for the rest of the 
variables (i.e. for fc > 2) when /3 ^ is 



fij \Cj 



N 



-iT^- (17) 

n=2(n#fc) " k 

Using ([12]) and the hmit of ([TT]) for fc = 1 as /3 ^ 0, the 
energy of the corresponding state can be found as 



N 



E 



N 



2Ne 



f^.ic.f + E 



(18) 



n=2 



E 



2£i 



2-f3 
P-Vk 



N 

E 

n=l{^k) 



Vn - Vk 



(25) 



In the limit where the single particle energies are degen- 
erate, (PT|) - (P5|) approach to those found in Ref. [9|. 

Next, we provide exact solutions of the equations of 
Bethe ansatz for iV = 1 and N — 2 pairs to illustrate 
the technique. Since the occupancy amplitudes are nor- 
malized as |cjj2 + |cj2p = 1 we can set cj-^ — sini) and 
= cos-d where < d < tt. For one pair, one has to 
solve the Bethe ansatz equation given in ([5]). This equa- 
tion has two distinct solutions which lead to two different 
energy eigenvalues when substituted in ([7]): 
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N=l 



51 j J sin^ -d - 



VLj^ cos^ d 



{Ej^ + Ej.,) ± -\/4A2 - 4 sin^ i9 - 17^, cos2 z9) A + (rj^-^ sin^ + fl^, cos^ ■d)' 



r 



(26) 



where we defined A = Sj^ — ej^ ■ 

For N = 2, one should solve the equations of Bethe 
ansatz (jlip for xi and X2- For two pairs, there are three 
distinct solutions of the equations of Bethe ansatz leading 
to three energy eigenvalues. Since the resulting analytical 
expressions for the energy eigenvalues are long, we find it 
more convenient to present them in an alternative form. 
Starting from the Bethe ansatz equations pT|) and the 



J 



energy it can be shown that for = 2 the three 

energy eigenvalues of the Hamiltonian are given by 



En= 



w 



cos(2i?) 



(27) 



where w is one of the three roots of the following cubic 
polynomial: 



2w^ + 2 cos 2z9 (3P ~2)w^ + cos^ 2z9 [4P (P - 1) + 8 A {^j, sin^ i} - cos^ i9) - 8 (A + sin^ i}) (A - cos^ i?) 
- sin^ 2i? {flj, + VLj^ )] w - 8 a cos^ 2d [Vtj^ {VLj^ - 1) cos'^ d - VLj^ (Q^, - 1) sin"* + AP] = 0. (28) 

I 



flj^ cos^ -d. Roots of 



Here we defined P 
this polynomial can be found analytically for any set of 
parameters. For example, in Figure [l] we plot the ex- 
act energy eigenvalues for two pairs found from (j27|) and 
([28|) for jl — 3/2 and j2 = 5/2. We show how the energy 
eigenvalues change with the occupation amplitudes for 
different values of A = Sj-^ — Sj^ ■ In the degenerate case 
(i.e. when A = 0), one of the energy eigenvalues is con- 
stant in agreement with (|16p . When we slowly deviate 
from the degenerate case, for example when A — ±0.2, 
this energy eigenvalue is still independent of the occupa- 
tion amplitudes for the most part. 

In this Letter, we presented the exact solution of the 
orbit dependent pairing problem between two nondegen- 



erate energy levels. We showed that this solution repro- 
duces the previously known results in the reduced pairing 
limit and the degenerate limit 0,BS[i3|- Since the solu- 
tion is exact, it is also valid away from these limits and it 
can be used to study the interplay between the one-body 
effects and the orbit dependent pairing effects in a realis- 
tic model. We also believe that the technique presented 
here is an important step towards the exact analytical 
solution of the most general pairing problem. 
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